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CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This application is related to co-pending application serial number (Attorney 

Docket No. 1391-268200) titled "Acoustic Logging Apparatus And Method." 

STATEMENT REGARDING FEDERALLY SPONSORED 
RESEARCH OR DEVELOPMENT 

[0002] Not applicable, 

BACKGROUND OF THE INVENTION 

Field of the Livention 

[0003] The present invention generally relates to acoustic well logging. More particularly, the 
present invention relates to determining slow and fast shear wave velocities and orientations in an 
anifiGtropic earth formation. 

Description of the Related Art 

[0004] It is a well known phenomena that certain earth formations exhibit a property called 
"anisotropy", wherein the velocity of acoustic waves polarized in one direction may be somewhat 
different than the velocity of acoustic waves polarized in a different direction within the same earth 
formation. See, for example, S. Crampin, A Review of the Effects of Anisotropic Layering on the 
Propagation of Seismic Waves, Geophys. J, R. Aslx. Soc, vol. 49, pp 9-27, 1977, incorporated 
herein by reference as if reproduced in foil below. Anisotropy may arise from intrinsic structural 
properties, such as grain alignment, crystallization, aligned fractures, or from unequal stresses 
within the formation. Anisotropy is particularly of interest in the measurement of the velocity of 
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shear/flexural waves propagating in the eartii formations. Shear or S waves are often called 
transverse waves because the particle motion is in a direction ^transverse", or perpendicular, to the 
direction that the wave is traveling. 

[0005] Acoustic waves travel fastest when the direction of particle motion (polarization direction) 
is aligned with the material's stiLffest direction. If the formation is anisotropic, meaning that there is 
one direction that is stiffer than another, then the component of particle motion aligned in the stiff 
direction travels fester than the wave component aligned in the other, more compliant direction in 
the same plane. A shear wave induced into an anisotropic formation splits into two components, 
one polarized along the formation's stiff (or fast) direction, and the other polarized along the 
formation's compliant (or slow) directioiL Generally, the orientation of these two polari2ations is 
substantially orthogonal (components which are at a 90° angle relative to each other). The fast 
wave is polarized along the direction parallel to the fracture strike and a slow wave in the direction 
perpendicular to it. 

[0006] Acoustic well logging techniques have been devised for determining the amount of 
anisotropy from the shear wave velocities (slowness), and the amojunt of anisotropy is generally 
defined as the difference between the velocities of the fast and the slow shear waves. One method 
of deterrnining fast and slow shear wave velocities and orientations uses an acoustic logging 
tool 100, as shown in Figure 1, to detect components of the acoustic signals at several levels of 
dipole receivers. See, for example U.S. Pat. No. 5,712,829 (hereinafter "the '829 patent") issued to 
Tang et al., incorporated herein by reference as if reproduced in Ml below. 
[0007] In the '829 patent, two dipole sources X and Y, 102, are oriented orthogonal to each other. 
Signals detected by the dipole receivers A 104, parallel to the X source, are referred to as XA 
signals when the X source is triggered. Similarly, signals detected by dipole receivers B 106, 
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parallel to the Y source 102, are referred to as YB signals when the Y source is triggered. Cross- 
component signals can also be detected by the perpendicular receivers when each soiarce is 
energized, and these signals are referred to as the XB and YA signals for the X and Y sources 
respectively. Thus, a total of four sets of signals are created for each dipole receiver pair for each 
set of firings of the sources X and Y. 

[0008] Each of the four sets of signals can be represented as a time series, each of which consists of 
a series of numbers indexed with respect to increasing tkne from the instant at which the respective 
source is energized. The abscissa value in each series of numbers represents amplitude of the 
received signal. It must be understood, however, that the signal received by any particular receiver, 
regardless of which transmitter was fired, contains information about both the fast and the slow 
waves. Stated otherwise, the signal received by any particular dipole receiver is a combination of 
the signal induced by the fast wave and the signal induced by the slow wave. Determining the 
slowness of flie fest and slow waves involves separating the fast and slow signals from tiie actual 
received signals. Various solutions to determine the fast and slow waveforms from the received 
signals incorporating both exist, for example, in U.S. Patent No. 4, 817, 061 issued to AJford et al., 
incorporated herein by reference as if reproduced in fiiU below. Once the fest and slow waveforms 
are decomposed from the composite received waveforms, prior art acoustic determinations are 
made as to Ihe slowness of each of the waves, hi particular, this slowness determination typically 
involves determioing a coherence/semblance of the decomposed waveforms. 
[00091 While semblance may create visually pleasing results, determinmg slowness in this matter is 
unsuitable for enror estimation. Consequently, an improved method to determine fast and slow 
shear wave velocity and oriaitation in an anisotropic formation is desired. 
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SUMMARY OF SOME OF THE PREFERBED EMBODIMENTS 
[0010] The preferred embodiments of the present invention comprise a method and apparatus for 
determining the slowness and orientation of the fast and slow shear waves ia an anisotropic earth 
formation. The apparati^ for making this determination preferably comprises two dipole 
transmitters, oriented substantially perpendicular to each other, mounted on a tool and designed 
for imparting acoustic energy into the surrounding formation. The tool further comprises a 
plurality of dipole receiver pairs, the receiver pairs spaced apart from each other and from the 
dipole transmitters. The dipole receivers in each dipole receiver pair are preferably oriented 
substantially perpendicular to each other. The preferred method of operation involves firing each 
dipole transmitter at each depth level sequentially, and obtaining a plurality of received 
composite signals with the dipole receivers, as the tool is slowly raised or lowered in the 
borehole. Each, of the received signals is a composite signal containing information about the 
fast and slow shear waves. Each receiver pair on the same elevation creates four received signals 
for each set of transmitter fixings. 

[0011] A plurality of transfer functions of the formation are assumed and a series of source 
waveforms or wavelets are estimated using the received waveforms and the assumed transfer 
functions. More particularly, the preferred embodiments assume a transfer function for the 
formation at issue, and then estimate, usmg each set of received signals, a series of source signals 
that created the received signals based on.the assumed transfer function. An objective function is 
created which is indicative of the similarity of the estimated source signals. Because the acttial 
source signals are preferably the same, a low value of the objective fonction indicates that the 
assumed formation transfer function was close to the actual formation transfer function. The 
source estimation preferably is repeated using multiple transfer fimctions (assumed strike angles 
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and slowness values). The values of the objective fimction calculated axe preferably plotted in a 
starting time verses slowness verses strike angle graph, with the strike angle being the ordinate, 
the slowness being the abscissa, and the starting time being the Z axis coordinate. Thus, for a 
series of assumed transfer functions, all at a particular single strike angle, a vertical plane of 
information is created. The process is repeated for a series of assumed strike angles ranging fi:om 
-90'" to +90'' (for a total of 180^), and at a plurality of slowness values within each assumed 
transfer function. From minima^ in the graph, the orientations of the fast and slow axis may be 
determined, the difference in slowness between the fast and slow waves maybe determined, and the 
error of the slowness calculation determined. 

[0012] The disclosed device comprises a combination of features and advantages which enable it to 
overcome the deficiencies of the prior art devices. The various characteristics described above, as 
well as other features, will be readily apparent to those skilled in the art upon reading the following 
detailed description, and by referring to the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0013] For a more detailed description of the preferred embodiments of the present invention, 
reference will now be made to the accompanying cirawings, wherein: 
[0014] Figure 1 shows a prior art acoustic logging tool; 

[0015] Figure 2 shows a wirehne acoustic Jogging tool of the preferred embodknent; 

[0016] Figure 3 shows an exemplary set of decomposed received signals; 

[0017] Figure 4 shows an exemplary plot indicating how objective function values are placed in the 

plots of the preferred embodiment; and 
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[0018] Figure 5 shows a black and white exemplary plot of an objective function against slowness 
and strike angle. 

NOTATION AND NOMENCLATURE 
[0019] Certain terms are used throughout the following description and claims to refer to 
particular system components. This document does not intend to distinguish between 
components that differ in name but not function. In the following discussion and in the claims, 
the terms "including" and "comprising' are used in an open-ended fashion, and thus should be 
interpreted to mean "includuig, but not limited to. . .". 

DETAILED DESCRIPTION OF THE PREFERRED ElMBODMENTS 
[0020] Figure 2 shows an acoustic logging device 10 constructed in accordance with the preferred 
embodiments. In particular, Figure 2 shows the tool 10 disposed within a fluid filled borehole 12 at 
some distance below the surface 14. The tool 10 is preferably suspended within the borehole by 
means of a multi-conductor armored cable 16. Thus, the tool 10 of the inrefenred embodiment is a 
wireline device. However, it must be understood that the principles described herein may likewise 
be practiced in a measuring-while-drilling (MWD) or logging-while-driUing (LWD) system. 
[0021] The tool 10 preferably comprises a set of dipole transmitters: a first dipole transmitter 20, 
and a second dipole transmitter 22. In the perspective view of Figure 2, only one face of each of the 
dipole transmitters 20, 22 may be seen. However, one of ordinary skill in the art understands that a 
complimentary face of each dipole transmitter 20 and 22 is present on a back surfece of the tool 1 0. 
The dipole transmitters each may be a single dipole transmitter extending through the center of the 
tool 10, or may be individual transmitters fired in such a way as to act in a dipole feshion. The 
transmitter 20 induces its xoustic energy along an axis, which for convenience of discussion is 
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labeled X in the Figure 2. Transmitter 22 preferably induces energy along its axis labeled Y in 
Figure 1, where the X and Y axes (and therefore traosmitters 20, 22) are preferably orthogonal. The 
orthogonal relationship of the transmitters 20, 22 need not necessarily be the case, but a deviation 
from an orthogonal relationship complicates the decomposition of the waveforms. 
[0022] Tool 10 also comprises a plurality of receiver pairs at elevations spaced apart from the 
transmitters 20, 22, In particular, the preferred embodiment comprises four pairs of dipole 
receivers. However, any number of receiver pairs may be used without departing from the spirit 
and scope of the invention. In Figure 2, the receivers are labeled 24A-D and 26A-D, Preferably, 
each set of dipole receivers at a particular elevation has one receiver whose axis is coplanar with the 
axis of transmitter 20 (in the X direction) and one receiver whose axis is coplanar with the axis of 
transmitter 22 (in the Y direction). For example, one set of dipole receivers could be receivers 24A 
and 26A. Thus, the dipole receivers whose axes are coplanar with the axis of transmitter 20 are the 
transmitters 24A-D, Likewise the dipole receivers whose axes are coplanar with the axis of 
transmitter 22 are receivers 26A-D. While it is preferred that the axes of the receivers be coplanar 
with the axes of one of the transmitters, this is not required. However, azimuthly rotating any of the 
receiver pairs complicates the trigonometric relationships and, therefore, the data processing. 
[00231 Idl broad terms, the processing method of the preferred embodiment comprises calculating 
or estimating source signals or source wavelets that created each set of received signals by 
assuming a transfer function of the formation. Estimating source wavelets can be described 
mathematically as follows: 

where S^si^ is the estimated source signal calculated for the ith set of receivers, XrF\ is the assumed 
transfer function of the formation in the source to receiver propagation, and i?/ is the decomposed 
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waveforms (described below) for the ith receiver set. Thus, for each set of received signals Ri, an 
estimate of the source signal S^st. created. The estimated sovarce signals are con^ared in some 
way to create an objective ftmctioii. Minimas of a graph of the objective fbnction are indicative of 
the angle of the anisotropy, and the slowness of the acoustic waves through the formation. Further, 
depending on the type objective function used, one or both of the value of the objection function at 
the minicnas, and the curvature of the of the objective function plot near the nainimas, are indicative 
of the error of the slowness determination. . 

[0024] Thus, a primary component of the source signal estimation is the assumed transfer function 
[TF]. The transfer function may be relatively simple, taking into account only the finite speed at 
which the acoustic signals propagate and the strike angle, or may be very complex, to include 
estimations of attenuation of the transmitted signal in the formation, paths of travel of the acoustic 
signals, the many different propagation modes within the formation i&g. compressional waves, 
sheer waves, Stonely waves), and if desired even the effects of the acoustic waves crossing 
boundaries between different layers of earth formations. For reasons of sknplicity of the 
calculation, the preferred estimated transfer functions take mto account only the propagation speed 
(slowness) of the acoustic energy in the formation and the strike angle of tiie anisotropy. 
[00251 As discussed in the Background section, anisotropic earth formations tend to break an 
induced shear wave into two components: one of those components traveling along the fester 
polarization direction, and the second conjponent traveling along the slower polarization direction, 
where those two directions are substantiaUy orthogonal. The relationship of the fast and slow 
polarizations within Ihe formation, however, rarely lines up with the orihogonal relationship of the 
dipole transmitters 20, 22. For convenience of the following discussion and mathematical 
formulas, a strike angle 9 is defined to be the angle between the X direction orientation (the axis of 
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dipole transmitter 20) and the faster of tlie two shear wave polarizations (see Figure 2). Further, it 
must be imderstood Uiat the shear wave of interest does not propagate in the X or Y direction, but 
instead propagates in the Z direction along the borehole wall. 

[0026] Operation of the tool 10 involves alternative fixings of the transmitters 20, 22. Each of the 
receivers 24A-D and 26A-D create received waveforms designated R, startii^ at the firing of a 
particular transmitter. Each of the received waveforms or signals has the following notation: 
R[rec«v«][sou«»]. Thus, for the firing of transmitter 20 in the X direction, and receipt by one of the 
receivers having an axis coplanar to the axis of transmitter 20 (receivers 24A-D), the time series 
received signal is designated as.Rxx. Likewise, the cross-component signal, the signal received by 
the dipole receiver whose axis is substantially perpendicular to the axis of the firing transmitter, is 
designated Rtx in this situation. In similar fashion, firing of the transmitter whose axis is oriented 
in the Y direction, transmitter 22, results in a plurality of received signals designated as Ryy for the 
axially parallel receivers, and Rxy for the cross-components. Thus, each transmitter firing creates 
two received signals, one for each receiver of the dipole receiver pair. It follows that for a set of 
dipole transmitter fffings, four signals are received at each receiver pair indicative of the acoustic 
signals props^ted through the formation. 

[0027] Each of the received signals in the case described above contains components of both the 
fest and slow shear waves ~ composite signals. That is, for example, an Rxx receiver signal 
contains information regarding both the fest and slow polarized signals. These composite signals 
may be decomposed into their fest and slow components using equations as follows: 

FP(0=cos^ (e) Rxx(£l + sin(e)cos(0) tRxy(0+R^(t)] +sin\e) Ryy (?) (2) 
SP(0 = sin' (9) Rxx it) - cos(e )sin(e) \R^t) + RYx(r)] + cos' (6) ^T({t) 0) 
sin(2e)[Rxx(0 - Ryy(0] - cos (28) [Rxy(0 +Ryx(0] = 0 <4) 
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where FP(t) is the fast polarization time series, SP(t) is the slow polarization time series, and 6 is 
the strike angle as defined above. The prior art technique for decomposing the multiple received 
composite signals involved determining the strike angle 0 by solving equation (4) above, and using 
that strike angle in equations (2) and (3) to decompose the composite signals into the fast and slow 
time series. The preferred embodiments of this invention take a different tact. 
[0028] A close inspection of equations (2) and (3) above for the fast and slow polarization time 
series respectively shows two very symmetric equations. Taking into account the trigonometric 
relationships: 

sme = cos(90^-e) (5) 
cos 9 = sin (90^ -G) (6) 
it may be recognized that either the fast polarization equation (2) or the slow polarization equation 
(3) may be used to obtain either the fast or slow polarization signals by appropriately adjusting the 
angle 0 used in the calculation. Stated otherwise, either the fast or slow polarization equations (2) 
or (3) may be used to decompose a received signal having both fast and slow components into 
individual components if the strike angle 0 is appropriately adjusted. 

[00291 Rather than using a single strike angle in both equations (2) and (3) above, in the 
preferred embodiments each assumed transfer function comprises a strike angle. A plurality of 
transfer functions are assumed over the course of the slowness determination, and thus a plurality 
of strike angles are used, preferably spanning possible strike angles from -90° to +90° (180°). 
For each assumed transfer function (and Ihus strike angle), the four received signals generated by 
a set of receivers at each elevation are decomposed using the following equation: 

DS (t) = cos' (0) Rxx(/) + sin(0)cos(0) [RxY(/)+R^(t)] W(e) Ryy (f) 0) 
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where DS (t) is simply the decomposed signal for the particular strike angle used. This process is 
preferably repeated for each set of received signals at each level for each assumed transfer 
fimction. Equation (7) is equation (2) above; however, equation (3) may be equivalently used if 
the assumed strike angle is appropriately adjusted. 

[0030] Figure 3 shows an exemplary set of four decomposed signals that in the preferred 
embodiment are created using equation (7) above for a particular transfer function (strike angle). 
In the exemplary set of decomposed signals, Rl could be the decomposed signal created using 
the strike angle from the assumed transfer function and the composite signals received by the set 
of receivers 24A, 26A. Likewise, decomposed signal R2 could be the decomposed signal created 
again using the strike angle from the assumed transfer function and the composite signals created 
by the set of receivers 24B, 26B. Notice how the amplitude of the decomposed signal of the set of 
receivers closest to the transmitters, decomposed signal Rl, is greater than the decomposed signals 
of the more distant receivers, for example R4. Note also how the waveforms shift out in time from 
the closest to the more distant receivers, which is indicative of the finite speed of the acoustic waves 
within the formation. 

[0031] For a particular starting time within the decomposed signals, for example starting time Ti in 
Figure 3, and for a first assumed transfer function having an assumed strike angle and slovraess, 
portions of each decomposed signal are identified as being related based on the transfer fimction. 
Rectangular time slice 50 of Figure 3 is representative of a slowness in an assumed transfer function 
(with the assumed strike angle used to create the deconq>osed signals exemplified in Figure 3). In 
particular, the slope of the rectangular time slice is indicative of the slowness of the assumed 
transfer fimction. Stated another way, the portions of the decomposed signals within the 
rectangular time slice 50 should correspond based on the assumed slowness of the fomiation of the 
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transfer function. The time width of the samples taken from each of the received sigaals is 
preferably at least as long as each of the source signals in a firing set In lliis way, an entire source 
waveform or source wavelet maybe estimated However, the time width of the samples taken firom 
the decomposed signals need not necessarily be this width, as shorter and longer times would be 
operational. 

[00321 Thus, the portions of the decomposed signals in the rectangular time slice 50 are each used 
to create an estimated source signal. These estimated source signals are compared to create an 
objective fimction that is indicative of their similarity. The process of assuming a transfer function, 
estibnating source wavelets based on decomposed signals and creating an objective function is 
repealed a plurality of times. The rectangular time slices 50 through 54 are exemplary of multiple 
assumed transfer functions used in association with starting time Ti (and the a strike angle used to 
create the decomposed signals). Estimating source wavelets in this fashion (includmg multiple 
assumed transfer functions) is preferably also repeated at multiple starting times within the 
decomposed signals, represented in Figure 3 as Ti, T2 . . . Tn- 

[0033] The value of the objective function created for each assumed transfer function and starting 
time is preferably plotted in a graph as a function of the starting time and the slowness of assumed 
transfer function. As shown in Figure 4, the starting time of the source signal estimations, Ti, 
T2 . . . Tn, are preferably the Z axis coordinate (assumed stiike angle is the ordinate or X axis (not 
shown in Figure 4), as is discussed more below) and the slowness is preferably the abscissa or Y 
axis. Thus, for example, the value of tiie objective function calculated using portions of 
decomposed signals Ri - R4 within rectangular time sUce 50 is plotted at point 60 in Figure 4. 
Likewise, the value of tiie objective fimction associated with the assumed slowness implied by 
rectangular time sUce 52 is plotted at point 62, and the value of the objective fimction associated 



54525.01/1391-27200 



-12- 



with the rectangular time slice 54 is shown at point 64. Thvs, for each of a pluraHty of starting 
times, and for each starting time a plurality of slowness values associated with assumed transfer 
functions, a graph of the objective function is created. In the preferred embodiments, the value of 
the objective function is indicated by a color, with cooler colors (blues, purples) showing more 
similarity, and hotter colors (reds, oranges) showing less similarity of the estimated source signals 
or source wavelets. The color schemes however are only exemplary and other schemes may be 
used without departing from the scope and spirit of the invention. Co-pending application Serial 

No. (Atty. Docket 1391-26800) titled, "Acoustic Logging Apparatus and Method", 

incorporated by reference herein as if reproduced in fuU below, describes related methods. 
[0034] Calculating the objective function of the a first embodiment preferably comprises 
comparing estimated source signals to determine a variance between them. More particularly, this 
slowness determination preferably comprises calculating an average of tiie estimated source signals 
within each time slice, and then calculating a variance against the average source signal. In more 
mathematical terms, for each assumed transfer fimction, a series of estimated source waveforms or 
signals S^st, are calculated using equation (1) above. From these estimated source signals, an 
average estimated source signal is preferably calculated as follows: 

where Sest is tibe average estimated source signal, JST is the number of decomposed received 
signals, is the source wavelet estimated for each decomposed received signal within the time 
slice, and t is time within the various time seri^. 
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[0035] The average estimated source signal is used to calculate a value representing the variance of 
the estimated source signals from the average estimated source signal. The variance is preferably 
calculated as follows: 

where ^ is the variance. In one embodhnent, the variance value is flie objective function plotted in 
the slowness versus starting time versus strike angle graph. Large values of the variance indicate 
that the ^sumed transfer fimction (assumed strike angle and/or assumed slowness) did not 
significantly match the actual formation properties. Likewise, small values of the variance indicate 
that the assumed transfer function closely matched the actual formation properties. Thus, the 
mmimas of the objective fimction in the plot described above indicate the slowness of the fest and 
slow polarized waves as well as the actual strike angle (this detennination further exemplified in the 
discussion with respect to Figure 5 below). The value of the variance objective fimction at the 
^inimas is indicative of the error of the determination of the acoustic velocity and strike angle. 
Relatedly, the curvature of the variance objective fimction plot at the minima is indicative of the 
error of the calculation. 

[00361 A second embodiment for calculating an objective fimction is based on determining a 
difference between each estimated source signal. In particular, and referring again to Figure 3, 
consider the portions of the decomposed signals within rectangular time slice 50 associating witii 
startmg time Ti. As discussed above, i^ing the assumed transfer function, an estimated source 
signal is created using the portions of the decomposed signal within the time slice 50. Differences 
or differentials are calculated between each estimated source signal, for example between the 
source signal estimated from a portion of the Rl signal and the source signal estimated from the 
portion of the R2 signal. This difference is preferably calculated between each succeeding receiver, 
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and the objective function in this embodiment is the sum of the square of each difference 
calculation. Much like plotting the variance objective function, this differential objective function 
is preferably plotted as a function of slowness , starting time and strike angle. However, the graph 
obtained using the differential slowness calculation has slower transitions from maximas to 
minimas which therefore makes determining the minimas (indicative of the actual slowness of the 
fast and slow polarizations) easier than in cases where the graph has relatively steep slopes between 
minima and maxima More mathematically, the objective function of this second embodiment is 
calculated as follows: 

i-l 

where ^ is the objective function, and N is the number of receivers. Much like using the variance 
as the objective function, this differential objective function is preferably plotted in a slowness 
versus starting time versus strike angle graph, with cooler colors representing less difference 
between received signals, and hotter colors representing greater differences. Known techniques 
may be used to determine minima of tiiese graphs, and the locations of the mmima are indicative of 
formation slowness and the strike angle. 

[0037] Neither of the two embodiments of calculating objective functions described above, 
variance and differential, are more preferred, as each have their own advantages. Plots of the 
variance objective value have steep slopes between maxima and minima, and thus tiie minima 
are better defined and the results more accurate. However, the differential objective value system 
has slower transitions from maxima to minima, makmg the determmation of the minima by 
computer program easier. Thus, either may be used m any particular circumstance, and it is 
possible that both may be used in a single system. 
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[0038] Figure 5 shows a black and white version of an exemplary slice (to create a two 
dimension graph) of a three dimensional graph of the objective function of the preferred 
embodiment. With reference to Figure 3, preferably slowness is the abscissa, starting time is the 
Z coordinate, and strike angle is the ordinate (extending out of the page). The plot of Figure 5 is 
a vertical slice of the three dimensional plot at a particular starting time, e.g. Tl, such what is 
produced is a two dimensional plot as a ftinction of slowness and strike angle. As can be seen 
from Figure 5, the slower of the two polarized shear waves has a slowness of approximately 140 
micro-seconds per foot ([is/ft), aad is polarized at approximately -45°, Likewise, the faster of the 
two polarized waves has a slowness of approximately 125 p,s/ft, and is polarized at 
approximately +45°. As discussed above, the two regions of lower objective function values are 
preferably plotted in color, with lower objective function values taking on cooler colors. 
Figure 5 however, being black and white, shows these only as darker shades of gray. By 
an^zing the three dimensional plot, using known techniques, to find minimas of the objective 
function, characteristics of the formation may be determined. These characteristics may 
comprise the fast and slow shear wave velocities, strike angle of the aoisotropy (as weU as the 
relative angle between the two waves), and an indication of the error in the calculation (the value 
of the objective function at the minimas). Another error indication is the calculated difference m 
angle between the fast and slow wave. Theoretically, the angle should be 90 degrees. A 
departure from the 90 degree theoretical value may be used as an estimate of the error of the 
strike angle calculation. 

[0039] Numerous variations and modifications will become apparent to those skilled in the art 
once the above disclosure is fully appreciated. For example, the disclosed method for 
determinmg shear wave velocity and orientation may be implemented using any number of 
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receiver levels and different receiver types for the acoustic logging tool. Indeed, even a single set 
of dipole receivers may be used relying on rotation of the tool to obtain additional composite 
signals for decomposition. Further, the source may be located at any arbitrary angle relative to 
the receivers. Moreover, processing of the data after collection at receivers can be performed 
downhole in real time with only the results being transferred uphole to a computer system for 
storage. Throughout this discussion, the various earth formation characteristics were discussed 
with reference to finding minimas of the objective function. However, one of ordinary skill in 
the art could easily invert the values used, thus making a determination a search for maximum 
values in the plot, and this would not deviate from the scope and spirit of the mvention. While 
assuming the transfer functions in the embodiments described involved also assuming a strike 
angle, it is possible that the transfer function need not include a strike angle estimation, and 
instead the composite signals could be decomposed for the range of possible strike angles 
in4gpendent of an assumed transfer fimction. Relatedly, it is possible to solve for the strike angle 
using equation (4) above and decompose the composite waveforms using that strike angle; and 
thereafter, estunate and apply transfer functions to the decomposed signals, thus also removmg 
the strike angle from the transfer function. It is intended that the following claims be mterpreted 
to embrace all such variations and modifications. 
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